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ABSTRACT 

We search for massive and compact galaxies (superdense galaxies, hereafter SDGs) at z = 0.03-0.11 in the 
Padova-Millennium Galaxy and Group Catalogue, a spectroscopically complete sample representative of the 
local Universe general field population. We find that compact galaxies with radii and mass densities comparable 
to high-z massive and passive galaxies represent 4.4% of all galaxies with stellar masses above 3 x 10 10 M Q , 
yielding a number density of 4.3 x 10~ 4 h 3 Mpc~ 3 . Most of them are SOs (70%) or ellipticals (23%), are red 
and have intermediate-to-old stellar populations, with a median luminosity- weighted age of 5.4 Gyr and a 
median mass-weighted age of 9.2 Gyr. Their velocity dispersions and dynamical masses are consistent with the 
small radii and high stellar mass estimates. Comparing with the WINGS sample of cluster galaxies at similar 
redshifts, the fraction of superdense galaxies is three times smaller in the field than in clusters, and cluster 
SDGs are on average 4 Gyr older than field SDGs. We confirm the existence of a universal trend of smaller 
radii for older luminosity- weighted ages at fixed galaxy mass. As a consequence, the median mass-size relation 
shifts towards smaller radii for galaxies with older stars, but the effect is much more pronounced in clusters 
than in the field. Our results show that, on top of the well known dependence of stellar age on galaxy mass, the 
luminosity-weighted age of galaxies depends on galaxy compactness at fixed mass, and, for a fixed mass and 
radius, on environment. This effect needs to be taken into account in order not to overestimate the evolution 
of galaxy sizes from high- to low-z. Our results and hierarchical simulations suggest that a significant fraction 
of the massive compact galaxies at high-z have evolved into compact galaxies in galaxy clusters today. When 
stellar age and environmental effects are taken into account, the average amount of size evolution of individual 
galaxies between high- and low-z is mild, a factor ~ 1 .6. 

Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: structure — galaxies: funda- 
mental parameters 



1. INTRODUCTION 

Recent observations have found that the sizes of high-z qui- 
escent and massive galaxies at z ~ 1 -2.5 are on average much 
smaller than those of high Sersic index (n > 2.5) galaxies in 
the local Universe (Daddi et al. 2005, Trujillo et al. 2006, Toft 
et al. 2007, Zirm et al. 2007, Buitrago et al. 2008, Cimatti et 
al. 2008, van Dokkum et al. 2008, van der Wei et al. 2008, 
Saracco et al. 2009, Cassata et al. 2011, Damjanov et al. 
201 1). Recent tantalizing observations from CANDELS may 
have identified the z > 2 compact, star-forming progenitors of 
the passive z = 1 -2 population (Barro et al. 2012). At a given 
stellar mass, z = 1 - 2 passive galaxies exhibit a wide range of 
sizes, from those with radii comparable to nearby galaxies to 
ultra-compact galaxies up to six times smaller (Saracco et al. 
2009, Mancini et al. 2010, Cassata et al. 2011). Searches for 
superdense galaxies in the Sloan at low-z have found very few 
such galaxies in the general field (Trujillo et al. 2009, Taylor 
et al. 2010). 

These results suggest a strong evolution in galaxy size be- 
tween z = 2 and z = 0. The galaxy size growth appears to 
be stronger for more massive galaxies (Ryan et al. 2012, 
Huertas-Company et al. 2012), and at least at z < 1 depen- 
dent on morphological type (Huertas-Company et al. 2012). 

Several mechanisms have been proposed to explain these 
observations. One such mechanism is expansion for mass 
loss, due to an active galactic nucleus (Fan et al. 2008) or stel- 
lar winds (Damjanov et al. 2009), though Ragone-Figueroa 
& Granato (201 1) have found that these are insufficient to ex- 



plain the observed evolution. Major mergers are another pos- 
sible process at work, but the major merger rate seems to be 
insufficient to explain the observed evolution. More frequent 
minor mergers are currently considered the most likely expla- 
nation for the evolution of galaxy sizes, since they reduce the 
effective stellar densities, strongly increasing the galaxy size 
while leaving the mass almost unaffected (Khochfar & Silk 
2006, Naab et al. 2009, Hopkins et al. 2009, Bezanson et al. 
2009, Hilz et al. 2012, Oser et al. 2012, see also Trujillo et 
al. 2011). However, several studies struggle to reconcile the 
minor merger scenario with the large amount of evolution re- 
quired when comparing at face value high- and low-z galaxy 
populations (Nipoti et al. 2009, 2012, Saracco et al. 2011, 
Newman et al. 2012, Oogi & Habe 2012). 

Two new aspects that have been recently highlighted may 
lead to a modification of the current paradigm of size evolu- 
tion. The first one has been the detection of a significant frac- 
tion of superdense galaxies with radii and masses comparable 
to high-z galaxies in nearby galaxy clusters (Valentinuzzi et 
al. 2010a, hereafter V10) and in clusters at intermediate red- 
shift (Valentinuzzi et al. 2010b). This may imply a strong 
environmental dependence of the distribution of massive and 
compact galaxies in the local Universe, but no study so far has 
carried out a detailed comparison between clusters and field 
at low-z. 

The second aspect concerns the effects of galaxy selection 
in the high-z samples. Several studies have found that smaller 
galaxies have older stellar populations. As a result, by select- 
ing galaxies that are already passive, high-z studies might be 
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preferentially selecting the most compact ones (Valentinuzzi 
et al. 2010, Saracco et al. 2009, 2011, Cassata et al. 2011, 
Szomoru et al. 2011). The existence of a trend of size with 
age at fixed mass is however still debated (e.g. cf. Valentin- 
uzzi et al. 2010 and Trujillo et al. 2012). 

In this paper we search for superdense galaxies and analyze 
the mass-size relation of galaxies in a new sample representa- 
tive of the general field population in the local Universe, the 
Padova-Millennium Galaxy and Group Catalogue (PM2GC, 
Calvi et al. 2011). This is the first non-Sloan-based work of 
this kind in the field at low-z, on a much smaller sky area than 
Sloan, but with a better quality imaging and higher spectro- 
scopic completeness. 

We investigate the incidence and properties of superdense 
galaxies as a function of environment, comparing the general 
field with clusters, and a range of finer environments includ- 
ing isolated galaxies and groups. We also focus on the de- 
pendence of the mass-size relation on galaxy ages and envi- 
ronment, attempting to quantify the importance of selection 
effects when inferring the size evolution of individual galax- 
ies. 

After presenting our dataset (§2), in §3 we show our re- 
sults for the fraction, number density and properties of su- 
perdense galaxies in the field, comparing with clusters in §4. 
The mass-size relation for different types of galaxies, differ- 
ent stellar population ages and as a function of environment 
is discussed in §5. We conclude by comparing with previous 
works (§6) and summarizing our results (§7). We release our 
surface photometry catalog describing it in the Appendix. 

Throughout this paper we will use the cosmology (Ho, Q m , 
= (70,0.3,0.7). All the masses in this paper have been 
scaled to a Kroupa (2001) IMF (0.1-100 M ). 

2. THE DATASET 

We draw our galaxy sample from the Padova-Millennium 
Galaxy and Group Catalogue (hereafter, PM2GC, Calvi et 
al. 2011) consisting of a spectroscopically complete sam- 
ple of galaxies at 0.03 < z < 0.11 brighter than M B < -18.7. 
This sample is sourced from the Millennium Galaxy Cata- 
logue (MGC, Liske et al. 2003, Driver et al. 2005), a B-band 
contiguous equatorial survey of ~ 3&deg 2 complemented by 
a 96% spectroscopically complete survey down to B=20. The 
image quality and the spectroscopic completeness are supe- 
rior to Sloan, and these qualities make it an interesting dataset 
to study galaxy sizes in a complete sample of galaxies at low- 
z. In particular, the MGC is based on INT data (2.5m tele- 
scope) obtained with a median seeing of 1.3" and at least 750 
sec of exposure, with a pixel scale of 0.333"/pixel, while the 
Sloan (again 2.5m telescope) has a median seeing of 1.5" in g 
(the closest band to our data), an exposure time of only 54. 1 
sec and 0.396"/pixel. As for completeness, 14% of all galax- 
ies and 25% of SDGs in our spectroscopic PM2GC sample do 
not have an SDSS spectrum. 

The whole PM2GC sample consists of 3210 galaxies and 
is representative of the general field population. Within the 
PM2GC, Calvi et al. (2011) identified 176 galaxy groups 
with at least three members brighter than Mb < — 18.7 using 
a Friends-of-Friends group-finding algorithm. A galaxy is a 
group member if its spectroscopic redshift lies within ±3cr 
from the median group redshift and if it is located within a 
projected distance of < 1 .57?2o(El from the group geometri- 
cal centre. Galaxies which have no neighbour or solely one 

1 /?2oo is defined as the radius of a sphere whose interior mean density is 



with a projected mutual distance < 0.5 h 'Mpc and a red- 
shift within 1500 km s" 1 are considered "single" galaxies or 
"binary-system" galaxies, respectively. 

The PM2GC sample is complete for masses M t > 1.8 x 
10 10 M Q (logM-t > 10.25), corresponding to the mass of the 
faintest and reddest galaxy (Mb = -18.7, B-V = 0.9) at our red- 
shift upper limit (z = 0.1), as described in Calvi et al. (2012). 
In this paper, however, we consider only those 1057 galaxies 
with > 3 x 1O 1O M (logM* > 10.48) to be consistent with 
the cluster work from Valentinuzzi et al. (2010). Surface pro- 
file fits were obtained for 94% of these galaxies as outlined 
below. Hereafter, we will use the mass-limited sample with 
radii estimates consisting of 995 galaxies at z=0. 03-0.1 1. 

2.1. Galaxy radii and Sersic indices 

Effective-radii, axial ratios and Sersic indexes are mea- 
sured on MGC B-band images with GASPHOT (Pignatelli & 
Fasano et al. 2006, Bindoni et al. in prep.), an automated tool 
which performs a simultaneous fit of the major and minor axis 
light growth curves with a 2D flattened Sersic-law, convolved 
by the appropriate, space-varying PSF In this way GASPHOT 
exploits the robustness of the ID fitting technique, keeping at 
the same time the capability (typical of the 2D approach) of 
dealing with PSF convolution in the innermost regions. 

GASPHOT has proved to be very efficient in recovering the 
best fitting parameters, and to give the appropriate weight to 
the external parts of the galaxies, where PSF effects are neg- 
ligible. Indeed we tested GASPHOT on more than 15,000 
simulated and real galaxies, obtaining robust upper limits 
for the errors of the global parameters of galaxies, even for 
non-Sersic profiles and blended objects (Pignatelli & Fasano 
2006). The AR e /R e error (median difference between re- 
trieved and input value) is -0.004±0.003 and -0.030±0.024 
for simulated n = 1 and n = 4 galaxies, respectively, and - 
0.057±0.101 for real galaxies. GASPHOT was also tested 
against the widely used tools GALFIT (Peng et al. 2002) and 
GIM2D (Marleau & Simard 1998). It was found (Pignatelli 
& Fasano 2006) that the performances of these tools are quite 
similar for large and regular simulated galaxies, while GAS- 
PHOT has proved to be more robust for real galaxies with 
some kind of irregularity, which is a crucial feature when deal- 
ing with blind surface photometry of huge galaxy samples. 

The GASPHOT output effective radius R e value is calcu- 
lated along the major-axis, and for the purposes of this paper 
is circularized with the usual formula: 

R (ciic) = R (majoi) . (1) 

where a and b are the major- and minor-axis of the best-fit 
model, respectively. 

Besides our own size measurements, there are four other in- 
dependent size estimates for PM2GC galaxies: GIM2D Ser- 
sic fits from the Millennium Galaxy Catalogue data (Allen 
et al. 2006), New York University Value Added Catalogue 
(NYUVAC) radii (Blanton et al. 2005), GIM2D Sersic fits 
from Sloan DR7 data (Simard et al. 2011), and DR7/SDSS 
de Vaucouleurs sizes accessed via the Catalog Archive Server 
(Abazajian et al. 2009, Thakar et al. 2008), the last three 
being based on Sloan g-band data. 

Figure [TJ shows the comparison between GASPHOT and 
GIM2D(MGC)+NYUVAC+SIMARD size estimates for the 

200 times the critical density of the Universe, and it approximates the cluster 
virial radius. 



Superdense galaxies at low redshift 





FIG. 1. — Top. Comparison between circularized GASPHOT radii and GIM2D/NYUVAC/SIMARD radii all in kpc. Bottom. Comparison between 
Vaucoulers DR7 radii (g band) and GASPHOT/NYUVAC radii as a function of GASPHOT radius and Sersic index. SDGs are red circles. 
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FIG. 2. — Comparison between the galaxy stellar masses used in this paper 
(PM2GC) and DR7 masses. Red galaxies are SDGs, black are galaxies in the 
mass-limited sample. 

618 galaxies in common among the four samples. The agree- 
ment is generally good, with a tendency for GASPHOT radii 
to be larger than the others. The median difference be- 
tween PM2GC radii and GIM2D/NYUVAC/SIMARD radii 
is 0.03±0.04/0.03±0.06/-0.01±0.04dex for all galaxies and 
0.01±0.08/-0.12±0.07/-0.06±0.11dex for superdense galax- 
ies (SDGs) as defined in §3. Therefore, for SDGs the Sloan- 
based estimates tend to be slightly larger than the MGC -based 
estimates, by a factor 1.3 and 1.1 for the NYUVAC and the 
Simard catalogues, respectively. The effect is clearly very 
small, especially when using GIM2D on Sloan data. 

Since there are previous low-z studies that have used de 
Vaucoulers DR7 radii (Taylor et al. 2010, Trujillo et al. 201 1), 
in Fig. [T] we also show the comparison between DR7 size es- 
timates obtained for a de Vaucoulers (n=4) fit and both GAS- 
PHOT and NYUVAC radii (764 in common among the three 
samples). The NYUVAC estimates have been obtained from 
the same data (Sloan g-band imaging) of the de Vaucoulers 
radii. As could be expected, the agreement is good for n ~ 4, 
but de Vaucoulers radii are much larger for lower n values, 
where most of the galaxies and, as we will see, most of our 
SDGs are found. 

2.2. Galaxy stellar masses 

Galaxy stellar mass estimates were derived by Calvi et al. 
(201 1) using the Bell & de Jong (2001) relation which corre- 
lates the stellar mass-to-light ratio with the optical colours of 
the integrated stellar population, using the B-band photometry 
taken from the MGC, and the rest-frame B-V color computed 
from the Sloan g — r color corrected for Galactic extinction 



(see Calvi et al. 201 1 for details). The typical uncertainty on 
stellar mass estimates is 0.2-0.3dex. 

Figure [2] shows that there is no systematic difference be- 
tween our mass estimates and DRc| ones, and the scatter is 
< O.ldex. The median mass difference (logM+(PM2GC)- 
logM+ (DR7)) is -0.01±0.07dex and 0.02±0.1dex for all and 
SDG galaxies, respectively. 

2.3. Morphologies 

To morphologically classify galaxies, we used MORPHOT, 
an automatic tool designed to reproduce as closely as pos- 
sible the visual classifications (Fasano et al. 2012). MOR- 
PHOT adds to the classical CAS (concentration/ asymme- 
try/clumpiness) parameters a set of additional indicators de- 
rived from digital imaging of galaxies and has been proved to 
give an uncertainty only slightly larger than the eyeball esti- 
mates (Fasano et al. 2011). MORPHOT was applied to the 
B-band MGC images to identify ellipticals, lenticulars (SOs) 
and later-type galaxies. These are the three broad morpholog- 
ical classes we will use in this paper (see Calvi et al. 2012 for 
details). 

It is well known that separating ellipticals from SOs could 
be sometimes challenging even for human experienced classi- 
fiers. The challenge is obviously enhanced when dealing with 
compact objects like SDGs. That being stated, the MOR- 
PHOT ability of recognizing visually classified SO galaxies 
on MGC-like imaging has been shown to be of the order of 
~ 80% (see Figure 11 in Fasano et al. 2012; see also Figures 
6 and 16 therein), which is comparable to the level of agree- 
ment reached by different experienced visual classifiers. 

2.4. Stellar population ages 

The MGC spectroscopic database of PM2GC galaxies con- 
sists of SDSS, 2dFGRS and MGCz spectra (Driver et al. 
2005), the latter taken with the same instrument and setup 
of the 2dFGRS. Choosing always the highest quality spec- 
trum available, we preferred to use SDSS spectra when pos- 
sible (858 cases, 86% of the sample), alternatively 2dF- 
GRS spectra when available (93 cases) and MGCz spectra 
in the remaining cases. SDSS spectra are flux-calibrated, 
while we performed a relative flux calibration for all 2dF 
spectra using the response curve provided by the 2dF team 
(http:www.mso. anu.edu. au2dFGRS) and a refined, second or- 
der calibration using a curve derived by the ratio between the 
2dF spectra and the corresponding SDSS spectra for objects 
with both spectra available, as described in Cava et al. (2009). 

To derive the stellar populations ages, we fitted the spectra 
with the spectro-photometric model fully described in Fritz et 
al. (2007). All the main spectro-photometric features (such 
as the continuum flux and shape, and the equivalent widths 
of emission and absorption lines) are reproduced by summing 
the theoretical spectra of Simple Stellar Population (SSP) of 
12 different ages (from 3xl0 6 to ~ 14xl0 9 years). 

Dust extinction is allowed to vary as a function of SSP age, 
in a screen uniformly distributed in front of the SSP stars. 
The Galactic extinction law follows the Cardelli et al. (1989) 
scheme, with Ry = 3.1. As explained in detail in Fritz et al. 
(2007), for the fit we use a fixed metallicity, exploring three 
values: Z=0.004, Z=0.02 and Z=0.05. The adopted star for- 
mation histories refer to the model with the metallicity value 
that provides the lowest \ 2 - 

2 http:www.mpa-garching.mpg. deSDSSDR7Datastellar mass.html 
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SSP spectra are built using Padova evolutionary tracks and 
the observed MILES spectral library (Sanchez-Blazquez et al. 
2006) for ages older than 10 9 yr, complemented by the Jacoby 
et al. (1984) library for young SSPs, and in the UV and in- 
frared by means of the Kurucz theoretical library. Nebular 
emission is included, modeled with values that are typical of 
Hll regions: this significantly affects spectra of SSPs younger 
than ~ 2xl0 7 years. 

From our spectral analysis, it is possible to derive an esti- 
mate of the average age of the stars in a galaxy. Following 
the definition of Cid-Fernandes et al. (2005), we compute the 
luminosity-weighted (LW) age by weighting the age of each 
SSP composing the integrated spectrum with its bolometric 
flux. This provides an estimate of the average age of the stars 
weighted by the light we actually observe. A mass-weighted 
(MW) age is computed in a similar way: each SSP age is 
weighted with its mass value. The mass-weighted age is the 
"true" average age of the galaxy's stars. By construction, the 
MW age is always older than the LW age. 

3. RESULTS: SUPERDENSE GALAXIES IN THE FIELD AT LOW 
REDSHIFT 

In the bottom panel of Figf3] we plot the circularized effec- 
tive radius R e as a function of stellar mass for PM2GC galax- 
ies with stellar masses M* > 10 1048 Mq. In the upper-panel, 
we plot the mean mass surface density inside P^: 



as usually defined by other authors (see, among others, 
Cimatti et al. 2008 and van der Wei et al. 2008). The plots 
include the median PM2GC mass-size relation for all (green, 
with lcr and 2a dashed lines), early-type (ellipticals + SOs, 
red line) and late-type (blue line) galaxies. The binned values 
of these mass-size relations are given in Table 1 . 

High-z literature data at 1 < z < 2.5 are plotted as large 
empty symbols: HUDF (Daddi et al. 2005, exagons), MU- 
NICS (Trujillo et al. 2006, circles), MUSYC (van Dokkum et 
al. 2008, squares), Saracco et al. (2009) (diamonds), GMASS 
(Cimatti et al. 2008, triangles), van der Wei et al. (2008) (pen- 
tagons), Damjanov et al. (2009) (reversed triangles), Mancini 
et al. (2010) (stars), Cassata et al. (201 1) (crosses 1 .2 < z < 2, 
starred crosses z > 2). All data have been converted to a 
Kroupa IMF (Kroupa 2001) and Maraston (2005) models as 
done in V10. The data, methods of analysis, galaxy mass 
ranges and selection criteria differ from one study to another, 
but all these works, with the exception of van der Wei et al. 
(2008) that have used a visual morphological classification, 
have selected galaxies to be already passively evolving with 
old stellar populations at the redshift they are observed, based 
on SED spectral fitting, absence of emission lines, line index 
age dating and other spectropho tome trie techniques. 

Following V10, we define SDGs trying to match as much 
as possible the position of most high-z SDG data in Fig. [3] 
applying the following criteria: 

3xlO lo M <M* <4xlO n M Q (3) 
£ 5 (>>3xl0 9 M Q kpc- 2 (4) 

SDG candidates complying these criteria have been indi- 
vidually inspected, to check their size and mass estimates. A 
merging system has been excluded from the analysis. When 
the GASPHOT and GIM2D estimates significantly disagreed, 



a galaxy -by-galaxy analysis performed with iraf was used to 
decide which radius was the most accurate. Most of these 
cases are galaxies with a nearby galaxy or star. There were 
three galaxies for which the DR7 mass estimate agreed with 
the mass obtained using SDSS CAS magnitudes but disagreed 
with the PM2GC mass. For these, the SDSS CAS mass was 
adopted. This left a total sample of 44 SDGs. 

The fraction of SDGs among PM2GC galaxies with masses 
3 x 10 10 < Mi < 4 x 10 n M© is 4.4 ±0.7%. The correspond- 
ing number density of SDGs within the PM2GC volume is 
4.3 x 10" 4 h 3 Mpc- 3 . 

We note that if we used only literature information the frac- 
tion of SDGs would be only slightly lower, but still compara- 
ble within the statistical error, 3.7±0.7% using GIM2D radii 
and DR7 masses and 3.3±0.7% for NYUVAC radii and DR7 
masses. 

We also note that the fraction of PM2GC galaxies and 
PM2GC SDGs with a size measurement from the NYUVAC 
are 85% and 73%, and the fractions with DR7 masses are 86% 
and 75%, respectively. This supports the notion that the SDSS 
sample might be indeed slightly biased against compact mas- 
sive galaxies (less complete for compact than non-compact 
galaxies), as discussed in detail by other authors (e.g. Taylor 
et al. 2010). 

It can be seen from Figf3]that our SDG definition corre- 
sponds to galaxies with sizes more than 2a smaller than the 
median mass-size relation. The region at radii below and 
densities above the solid straight line is the "zone of avoid- 
ance" where galaxies are rare in the ATLAS3D sample of 260 
nearby early-type galaxies, a complete and statistically rep- 
resentative sample of the nearby early-type population (Cap- 
pellari et al. 2011). Our SDGs are truly exceptionally dense 
objects, and outliers in the local mass-radius and mass-density 
relations. 

In Fig. [4] we compare the distributions of galaxy mass, ra- 
dius and density values of the PM2GC SDG sample with 
those of the high-z SDG data, together with the cluster low- 
z WINGS sample that will be discussed later. We do this 
to check that in the local Universe we are selecting galaxies 
with a "degree of compactness" that is comparable to that of 
high-z galaxies. For masses, we plot the simple distributions, 
while for radii and densities we plot the offset with respect 
to the high-z median mass-radius and mass-density relations 
obtained using all the high-z samples considered in this paper. 

The PM2GC, WINGS and high-z distributions largely over- 
lap (cf. Fig. 7 in Taylor et al. 2010). However, the low-z and 
high-z distributions do not coincide, and they are not expected 
to, simply because the low-mass datasets are mass-limited 
galaxy samples, while the high-z histograms are the results of 
a collection of different surveys, with different mass ranges, 
selection criteria etc. As a result, the mass distributions at 
high- and low-z differ, with the high-z histogram being more 
populated at high masses than the local samples, and the con- 
sequent differences are visible also in the radii and density 
panels. Nevertheless, the large overlap of the ALogR e and 
ALogJ^so distributions of low-z and high-z SDGs shows that 
we are indeed comparing local and distant samples of galaxies 
with a large overlap in degree of compactness and density. 

We note however that z > 2 datapoints (yellow histograms 
in Fig. [4} have more prominent low-radii and high-densities 
tails with respect to the z < 2 distributions. This is mainly 
due to the van Dokkum et al. (2008) sample, while Cassata et 
al. (201 1) datapoints resemble the z < 2 distribution. It is un- 
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TABLE 1 

PM2GC MASS-RADIUS RELATIONS. THE VALUES ARE LOGARITHM OF 
THE MEDIAN ESTIMATES, ERRORS ARE THE LOWER AND UPPER 
QUARTILES OF THE DISTRIBUTIONS. 

clear whether the fact that z > 2 galaxies seem to be extreme 
in their compactness is due to rapid evolution between z = 2 
and 1 (Cimatti et al. 2012) or selection effects (Saracco et al. 
2010), but observationally the difference with z appears quite 
strikingly and is worth further investigation. 

3.1. The properties of superdense field galaxies at low-z 

23±7% of our SDGs are ellipticals, 70±13% are SOs and 
only 7±4% are later-type galaxies. Figure [5] shows some ex- 
ample postage stamps of PM2GC elliptical, SO and late-type 
SDGs. In the whole mass-limited sample, 29% are ellipticals, 
30% are SOs and 41% are late-typesQ 

In Fig. [6] we present the distributions of the most relevant 
quantities describing our SDG sample: axial ratio, (B-V) rest 
frame color, Sersic index, radius, LW age, MW age, local 
galaxy density, and stellar massQ For comparison, overplot- 
ted are also the distributions for non-SDG galaxies. Table 2 
summarizes the mean values for SDGs. 

Most of the SDGs lie on the color-magnitude red sequence. 
Their color distribution is similar to that of the whole mass- 
limited sample, showing that being red is a characteristic of 
galaxies in this mass range regardless whether they are com- 
pact or not. 

Our SDGs show the tendency to be flattened on average, 
with a < b/a >= 0.48, as might be expected given the high 

3 Comparing with the Sloan Bayesian automated morphological classifi- 
cations of Huertas-Company et al. (2011), our morphological classifications 
agree with theirs (Prtype > 0.45, our SOs versus their SO+early-spirals) in 
79% of the cases for the 33 SDGs in common. 

4 Here the galaxy local density has been computed from the circular area 
that, in projection on the sky, encloses the 5 nearest projected neighbours 
brighter than My < -19.85 within ±1000kms~', minimizing survey edge 
effects using the SDSS and 2dFGRS in the PM2GC surrounding regions (see 
Vulcani et al. 2012 for details). 



PM2GC WINGS B-band 



Quantity 


Value 


RMS error 


Value 


RMS error 


# SDGs 


44 


7 


51 


7 


(Re) 


1.45 


0.26 


1.57 


0.34 


(n) 


2.8 


0.6 


3.1 


0.8 


(b/a) 


0.48 


0.13 


0.65 


0.16 


<M,>M 


6.0X10 10 


1.9xl0 10 


9.1xl0 10 


3.6xl0 10 


(Vabs) 


-20.87 


0.42 


-20.68 


0.38 


(Lw-age) 


5.45 


1.87 


9.64 


2.10 


(Mw-age) 


9.25 


1.08 


11.95 


1.39 


Ellipticals frac. 


22.7% 


7.2% 


29.1% 


7.8% 


SOs frac. 


70.5% 


12.7% 


62.0% 


10.7% 


Late-type frac. 


6.8% 


3.9% 


8.8% 


4.4% 



TABLE 2 

Characteristic numbers of PM2GC SDGs, compared to 
WINGS. Errors are derived from Poissonian statistics for 
counts and fractions, and ARE RMS FOR OTHER quantities. 



number of SOs0 Similarly, most of our SDGs have Sersic 
indices lower than 4, with < n >= 2.8 and a distribution inter- 
mediate between disky and spheroidal galaxies, as also high- 
z compact quiescent galaxies do (Wuyts et al. 2010, van der 
Wei etal. 2011). 

Most SDG radii are between 1 and 2 kpc, with a mean 
< R e >= 1.4. Stellar masses can be as high as 10 1L6 , but 
most of them are below 10 11 with a mean 10 10 ' 8 (6 x 10 10 M Q ). 
SDGs are intrinsically rather bright, with a mean absolute V 
magnitude of -20.9. 

PM2GC SDGs have an intermediate median LW age 
(5.4±1.5Gyr) and a old median MW age (9.2±0.8 Gyr), 
showing that the bulk of the mass has an age typically 4 Gyr 
older than the luminosity-weighted age. The great majority 
(93%) of SDGs galaxies have LW age > 1.5 Gyr. We have 
also inspected their spectra to search for the presence of emis- 
sion lines, and found that only about 10% of the SDGs have 
a significant emission (> 2 - 3 A) in any of the strongest lines 
([Oil] or Ha), never stronger than 15 A in [Oil]. 

The issue of the variation of galaxy sizes with their stellar 
population ages is particularly relevant when investigating the 
size evolution of galaxies. In addition to the well-known in- 
crease of mean stellar age with galaxy mass, some works, and 
in particular the cluster low-z work of V 1 0, have found that, at 
fixed mass, smaller (more compact) galaxies are older. This is 
confirmed by several works both at high- and low-z (Saracco 
et al. 2009, Shankar & Bernardi 2009, Kriek et al. 2009, 
Williams et al. 2010, Cappellari et al. 2012b), although for 
example Trujillo et al. (201 1) find no such a trend. Since most 
of the high-z studies select massive galaxies that are passive, 
whose stellar populations are already quite old at high-z, they 
would tend to select the smallest, i.e. more compact galaxies, 
with a strong selection bias (V10, Saracco et al. 201 1). 

To quantify to what extent, on average, PM2GC SDGs are 
older than non-compact galaxies of similar masses, we used 
the Monte carlo technique to extract 1000 random samples 
of "normal" galaxies with the same mass distribution of the 
SDGs. In Fig. [7] we plot the distributions of the median LW 
and MW ages of these 1000 samples. The Monte Carlo simu- 
lation shows that compact galaxies on average are ~ 0.7 and 
~ 0.2 Gyr older in LW and MW age, respectively, than nor- 
mal galaxies of the same mass (the median LW and MW ages 

5 The average b/a of high-z galaxies, found using Daddi et al. (2005), van 
Dokkum et al. (2008) and Damjanov et al. (2009), is 0.63±0.23. 
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FIG. 3. — The circularized effective radius R e and the mass-density inside R e as a function of stellar mass for all PM2GC galaxies with M» > 10 la48 M Q . The 
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line is delimiting the ATLAS3D zone of avoidance (see text). The green (with dotted la and 2<r lines), red and blue lines are the PM2GC median relations for 
all galaxies and for early and late-type galaxies, respectively. Open symbols are galaxies from high-z studies, see text for references. 



for normal galaxies are 4.7±1.8 and 9.0±0.8 Gyr). Again, 
this is suggesting that in some way age is related to compact- 
ness, in addition than mass, although the effect is modest for 
the LW-age and barely appreciable for the MW age. 

We will come back to this point and compare with LW and 
MW ages of cluster SDGs in §4, discussing the relevance of 
this effect for evolutionary studies. 

3.2. Velocity dispersions and dynamical masses 

As a further check of our SDGs mass and density estimates, 
we use the galaxy central velocity dispersions available for a 
subset of 570 galaxies in our mass-limited sample from the 
Sloan DR7 database. We compute the dynamical mass of 
all early-type (ellipticals + SOs) galaxies in our mass-limited 
sample following the formulation by Bertin et al. (2002): 

M dyn = (K(n)a 2 R e )/G (5) 

where 

/sT(«) = [73.32/(10.465 + (n-0.95) 2 )] + 0.954 (6) 

using the GASPHOT Sersic index n and circularized ra- 
dius R e and the DR7 velocity dispersion a corrected to R e 
as suggested in Cappellari et al. (2006). Using a fixed K(n)=5 
as recommended by Cappellari et al. (2006,2012) does not 
change the results. 



The comparison of virial and stellar mass estimates is 
shown in Fig. [8] For the great majority of galaxies, as well as 
the great majority of SDGs, the dynamical mass is higher than 
the stellar mass by a factor that is very similar to values found 
in the literature (Taylor et al. 2010b). The median difference 
logM+- logMjvn is -0.24dex. There are very few cases with 
stellar mass estimates higher than dynamical estimates, com- 
patible within the errors. The comparison between dynamical 
and stellar masses confirms the validity of the mass estimates 
we have used so far and confirms there is no problematic issue 
regarding our mass estimates. 

As a further, final test we check whether the velocity dis- 
persions of SDGs are higher than non-SDG galaxies of sim- 
ilar mass, as expected if they are indeed more compact than 
the normal population. Figure [9] shows that our SDGs (red 
symbols) have significantly higher velocity dispersions than 
the average non-compact galaxies of similar stellar mass. The 
median relations for SDGs and non-SDG galaxies are also 
shown. For SDG galaxies, the median offset of logR e from 
the mass-radius relation is -0.4dex, while the median Aloga 
is 0.1 dex. These values are very similar to the correspond- 
ing ones found by Taylor et al. (2010) for their SDSS sample 
(-0.38 and 0.12dex, respectively). The dynamical mass esti- 
mates and the velocity dispersions of our SDGs are therefore 
consistent with the fact that these are massive and compact 
objects. 
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FIG. 4. — Distribution of stellar mass, offset of effective radius from the high-z median mass-size relation and offset of density £50 from the high-z median 
mass-density relation. Symbols: PM2GC SDGs (black histogram), WINGS SDGs (blue shaded histogram) and high-z samples at z < 2 (red histogram) and > 2 
(yellow). Only high-z galaxies fulfilling SDG criteria (with 3xlO lo M Q < < 4xlO u M and S 50 > 3xl0 9 M Q kpc- 2 ) are included in this plot. Note that the 
median mass-size and mass-density relations used here include all galaxies at high-z, and this is the reason why the red+yellow histograms are not centered on 0. 



4. ENVIRONMENT 

As explained in §2, within the PM2GC survey, Calvi et al. 
(2011) have identified different environments: single galax- 
ies, binary systems, and groups with a wide range of velocity 
dispersions. These environments span a range of dark mat- 
ter halo masses, of the order of 10 12 5 to 10 15 M Q . We can 
therefore quantify the occurrence of SDGs in different envi- 
ronments. 

We find that the SDG fraction is always around 4%±1%, 
hence indistinguishable, in groups, binaries and single galax- 
ies. However, the fraction rises to 14±7% in PM2GC sys- 
tems with a velocity dispersion > 500 km s -1 . Splitting groups 
further into < 300 kms" 1 and SOO-SOOkms" 1 we do not de- 
tect a difference in SDG fraction, but uncertainties are large 
and there might be a more gradual transition in SDG fraction 
with environment which we do not detect due to small number 
statistics. 

In constrast with high-z results finding a trend of sizes with 
local galaxy overdensity (Cooper et al. 2012), there is instead 
no preference for SDGs to be in high density regions within 
the PM2GC, as shown in Fig.|6] 

4.1. Field versus cluster superdense galaxies 

SDG galaxies in nearby clusters were studied by the 
WINGS survey (V10). WINGS is a photometric and spectro- 
scopic multiwavelength survey of 77 X-ray selected galaxy 



clusters at z = 0.04-0.07 with a cluster velocity dispersion 
typically in the range 500-1200 kms" 1 and L x = 0.2-5 x 
lO^ergs" 1 (Fasano et al. 2006, Cava et al. 2009). 

Effective radii and galaxy morphologies were obtained 
from V-band images with the same tools used here, GAS- 
PHOT and MORPHOT, respectively (Pignatelli & Fasano 
2006, Fasano et al. 2012). Stellar galaxy masses were esti- 
mated both using the method adopted in this paper and with 
the spectrophotometric model used to derive stellar popula- 
tion ages, finding excellent agreement between the two meth- 
ods (Fritz et al. 201 1, Vulcani et al. 201 1). LW and MW ages 
were derived from optical spectra with the same method and 
spectrophometric model we used for the PM2GC (Fritz et al. 
2007,2011). 

The WINGS and PM2GC surveys, together, explore the 
whole range of environmental conditions, from massive clus- 
ters to groups, binary and isolated galaxies. All the tools 
and methods of analysis are the same for WINGS and the 
PM2GC, allowing a homogeneous comparison between mas- 
sive clusters and general field in the nearby Universe. 

For the purposes of this paper, the only relevant difference 
between the two surveys is that the size measurements in V10 
were performed on V-band WINGS images, instead of ri- 
band as in the PM2GC0 To obviate this problem, we have 

6 The consistency between V and B-band WINGS morphological classifi- 
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FIG. 5. — Example postage stamps of PM2GC elliptical, SO and late-type 
SDGs. 

employed GASPHOT to obtain effective radii from WINGS 
B-band imaging (Bindoni et al. in prep.), which was taken 
with the same instrument (WFCAM/INT) and similar depth 
of the MGC images we used in this paper. 

It turned out that the difference between the two filters is 
small, as shown in Fig. [TO] but should not be neglected. The 
median difference between the WINGS B-band R e and the V- 
band R e is 0.20 ±0.19 kpc for all galaxies, equivalent to a 
median AR e /R e ~ 0.1, and 0.29 ±0.30 kpc for SDGs with 
median AR e /R e ~ 0.17. GASPHOT was also run on WINGS 
K-band images (Valentinuzzi et al. 2009, Bindoni et al. in 
prep.) and we find a small but clear trend of galaxy sizes with 
wavelength as radii increase going from K to V to B. In the 
following, we only use B-band WINGS results. 

This small, but systematic offset has an effect on the SDG 
fraction. Adopting the usual criteria to define SDGs (equa- 
tions 3 and 4), B-band SDGs in WINGS are 11.8±1.7% of 
all galaxies more massive than our adopted mass limit. This 
fraction is about half of that found in the V band for the same 
limit (22±2%, V10). The B-band SDGs are essentially a sub- 
sample contained in the V-band SDG sample. 

The WINGS SDG proportion is consistent with the rough 
estimate obtained for > 500km s" 1 systems from the PM2GC 

cations is discussed in Calvi et al. 2012. 



Number density 
h 3 Mpc' 3 



General field SDGs (PM2GC) 

PM2GC SDGs quiescent today (LWage > 1.5Gyr) 

PM2GC SDGs quiescent z=1.0 (LWage > 9Gyr) 

WINGS "cluster" B-band SDGs 

WINGS B-SDGs quiescent z=1.0 (LWage> 9Gyr) 



4.3 ± 0.06 10 -4 
4.0 ±0.06 1CT* 
5.9 ±0.6 10~ 5 

> 1.7 ±0.09 10~ 5 

> 1.0 ±0.06 10 -5 



TABLE 3 

Number densities of SDGs and quiescent galaxies. Errors 
are derived from poissonian statistics. 



(14 ±7%, §3.1). SDGs are therefore proportionally much 
more numerous in clusters than in other environments, rep- 
resenting a significant fraction of the cluster massive galaxy 
population but a factor of 3 smaller fraction in the field. 

Let us now compare the number density of SDGs in the 
field and in clusters, summarized in Table 3. As said in 
§3.0, the number density of SDGs in the PM2GC is 4.3 x 
10" 4 h 3 Mpc" 3 . Considering only the oldest SDGs, those that 
have been quiescent at least since z=1.0 (=have stopped form- 
ing stars at least 1 Gyr before z=l, corresponding ~ to a LW 
age > 9 Gyr), the number density in the PM2GC is almost an 
order of magnitude smaller, 5.9 x 10~ 5 h 3 Mpc" 3 Q 

The WINGS B-band SDGs account for a number density 
of at least 1.7 x KrVMpcT 3 at z = 0.04-0.07. This num- 
ber has been computed by multiplying the average SDG num- 
ber per cluster among the 21 clusters considered in this study, 
corrected for spectroscopic incompleteness, by the total num- 
ber of clusters in the WINGS survey, and dividing by the 
whole volume of the Universe between z = 0.04 and z = 0.07. 
Considering the whole volume and not just the cluster vol- 
ume it is like assuming there are no SDGs outside of clusters, 
obtaining a hard lower limit on the SDG number density at 
these redshifts. We have not taken into account the fact that 
WINGS images are sampling the cluster populations as far as 
~ 0.6 x 7?2oo for all clusters, while they reach the virial radius 
for a small subset of clusters. Hence, the numbers above in- 
clude only galaxies in a fraction of the cluster volume, and the 
WINGS number density estimates are hard lower limits to the 
total SDG density accounted for by clusters. 

WINGS SDGs are generally very old, with a median LW 
age of 9.5Gyr, and the number density of LW age> 9 Gyr 
WINGS SDGs is at least (for the reasons given above) 1 .0 x 
10" 5 h 3 Mpc- 3 . 

Although we cannot obtain a precise estimate of the fraction 
of low-z SDGs that reside in clusters, due to the reasons out- 
lined above, we can conclude that, among old SDGs in the lo- 
cal Universe, a significant fraction (> 17%= 1.0 x 10~ 5 /5.9x 
10~ 5 ) must be found in clusters today. 

4.2. Cluster SDG properties 

Table 2 presents the average properties of WINGS B-band 
SDGs0 They have been computed weighting the values for 
spectroscopic incompleteness, as in V10. 

7 A direct comparison with published values of the number density of high- 
z compact galaxies is unfortunately unfeasible, because compactness criteria 
do not match and/or the largest high-z samples are affected by the z = 1.6 
overdensity in the GOODS-S field (Cassata et al. 2011, Barro et al. 2012). 
A careful comparison of the low-z and high-z number densities will be very 
valuable to assess what fraction of compact galaxies have survived till today 
(Renzini 2012). 

8 We note that these properties are similar to those of V-band SDGs given 
in Table 3 of V10. The distributions of SDG quantities are also similar to 
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FIG. 6. — Distribution of different quantities of interest for PM2GC SDGs (red solid histogram) and non-SDGs (black dotted): axial ratio, Sersic index, 
circularized effective radius, LW age, MW age, local density and stellar mass. The dashed black histogram in the Sersic index panel is the n distribution for the 
absolute magnitude-limited sample (see text). 
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FIG. 7. — Distributions of median luminosity-weighted ages (left blue his- 
togram) and mass-weighted ages (right red histogram) of 1000 random ex- 
tractions of non-superdense early-type galaxies with the same mass distribu- 
tion of SDGs. The vertical lines are the median ages of SDGs (solid lines) 
and of non-superdense galaxies (dashed lines) with their intrinsic mass dis- 
tribution. 
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FIG. 8. — The dynamical versus stellar mass relation for all PM2GC early- 
type (Es+SOs) galaxies (small black dots) and for SDGs (large coloured cir- 
cles - red for ellipticals and green for SOs). The typical errorbar is shown in 
the bottom right corner. The distribution of logarithmic difference between 
the two masses is shown in the inset for all galaxies and for SDGs (empty 
black and shaded red histograms, respectively). 



Fig. 6 in VI 0, and the latter can be compared with Fig. 6 in this paper. While 
the SDG fraction depends on the band used for their selection, the SDG prop- 
erties don't. 



FIG. 9. — Velocity dispersion versus stellar mass relation for non-SDG 
early-type galaxies (ellipticals + SOs, black) and superdense early-type galax- 
ies (red) for which a velocity dispersion measurement is available from Sloan. 
The black line is the median relation for non-SDGs, red for SDGs only. Ve- 
locity dispersions have all been corrected to R e as for Fig. [8] 

Their mean effective radius, Sersic index, axial ratio and 
absolute V luminosity are rather similar to the mean values of 
PM2GC SDGs. As in the field, the majority of cluster SDGs 
are SOs (62%) or ellipticals (29%). There is a tendency for 
cluster SDGs to be on average more massive than field SDGs 
(9 versus 6 x 10 10 Mq), but the most striking difference is in 
the stellar population ages. 

Cluster SDGs are on average 4 Gyr older than field SDGs 
in LW age, and about 3 Gyr older in MW age. The age dis- 
tributions for B-band WINGS SDGs are not dissimilar from 
those given in Fig. 6 of V10 for V-band selected SDGs and 
can be compared with Fig. 6 in this paper. Only 3 out of 51 
cluster SDGs have a LW age younger than 6 Gyr, while about 
half of the field SDGs are this young. 

The relation between galaxy mass, size and stellar age as a 
function of environment will be presented in detail in the next 
section, where we draw our conclusions regarding the evolu- 
tion of high-z massive compact galaxies and their descendants 
in the local Universe. 

5. THE MASS-SIZE RELATION: DESCENDANTS AND 
PROGENITORS 

In the previous section we have seen that the incidence of 
SDGs and the ages of their stellar populations are different in 
clusters and in the general field. In this section we wish to ad- 
dress a more general question: how does the relation between 
galaxy mass, size and stellar age depend on environment, and 
what can this tell us about the local Universe descendants of 
the high-z massive and compact galaxies. 

Let us first consider the mass-size relation. Figure QT| com- 
pares the median mass-size relations obtained in the PM2GC 
and in WINGS (B-band) for all galaxies and for galaxies of 
different types. We consider separately early-type (T-type< 0) 
versus late-type galaxies (T-type> 0), and galaxies with Ser- 
sic n > 2.5 versus n < 2.5. The binned values are listed in 
Tables 1 and 4. Figure [TTI includes also a comparison with 
the relations given by Shen et al. (2003) (« > 2.5 and n < 2.5) 
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FIG. 10. — Differences in circularized R e estimates in the B and V bands 
for WINGS galaxies. Red circles are V-band SDGs. The errorbar in the right 
upper corner represents the median relative error in the B-band AR e /R e =0.1. 



that are based on Sloan r-band imaging. The following con- 
clusions can be drawn from this figure: 

a) There is an offset between the mass-size relation in the 
general field and in clusters. Galaxies of a given mass in the 
field are bigger than in clusters typically by about la. This 
effect is visible for all galaxies, and also for each type of 
galaxies when dividing them by morphological type or Ser- 
sic index, being negligible only for late-type massive galax- 
ies. A KS-test can rule out the hypothesis of similar size dis- 
tribution for cluster and field galaxies of a given mass with 
probabilities 2 x 10" 8 at 10.5 < logM* < 10.8, 1.7 x 10" 8 
at 10.8 < logM* < 11.1, 1 x 10" 4 at 11.1 < logM+ < 11.4. 
Considering only early-type galaxies, the probabilities are 
6 x 10~ 3 , 9 x 10~ 5 and 3 x 10~ 5 in the three mass bins, respec- 
tively. We note that Huertas-Company et al. (2012b), instead, 
do not find a change in the median mass-size relation of early- 
type galaxies in the SDSS when comparing the general field 
and clusters. 

This is all based on B-band imaging, therefore it cannot 
arise from color gradient effects. We note that the WINGS B- 
band relation is already shifted to larger radii compared to the 
V-band relation for the same galaxies due to the wavelength 
dependence of effective radius shown in Fig.flOl 

b) The Shen et al. relation is shown here only to allow a 
comparison with the several works that have used it. When 
dividing by Sersic index, the PM2GC relation for n > 2.5 is 
shifted to larger radii compared to Shen's by about la, while 
the median relation is rather consistent for n < 2.5. We re- 
mind that an offset to larger radii is expected when using bluer 
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TABLE 4 

WINGS MASS-RADIUS RELATIONS IN THE B BAND. THE VALUES ARE 
LOGARITHM OF THE MEDIAN ESTIMATES, ERRORS ARE THE LOWER AND 
UPPER QUARTILES OF THE DISTRIBUTIONS. 



bands. Moreover, the Shen et al. work is known to be affected 
by systematic errors due to background oversubtraction (Guo 
et al. 2009), a bug in the NYUVAC size estimates used at the 
time that has been fixed in later NYUVAC releases (Blanton 
et al. 2005) and by the SDSS incompleteness at low redshifts 
(Taylor etal. 2010). 

c) Dividing galaxies by morphological types and by Ser- 
sic index is of course not fully equivalent. Differences can 
be observed in the mass-size relations derived for late-type 
galaxies (Tables 1 and 4). This is to be expected, given that 
dividing galaxies by Sersic index corresponds to include in 
the n < 2.5 group some early-type galaxies, that are gener- 
ally smaller than late-types and have the effect to decrease the 
average size value at a given mass. 

Let us now consider the effects of the stellar population age 
on the mass-size relation. In §3.2 we have seen that the me- 
dian LW age of SDGs is 0.7 Gyr older than the LW age of a 
sample of non-compact galaxies with the same mass distribu- 
tion of SDGs. In constrast, V10 found cluster SDGs to be on 
average 1.5 Gyr older than non-compact of the same mass. 

Figure Q~2] presents the median mass-size relation in the 
PM2GC (left) and B-band WINGS (right) for galaxies of dif- 
ferent LW ages. For the PM2GC we show galaxies younger 
than 4 Gyr (blue) and older than 8 Gyr (red), that represent 
the youngest and oldest wings of the age distribution. For 
WINGS, where the age distribution is skewed to much older 
ages, we can plot galaxies younger than 4 or 6 Gyr (blue), and 
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FIG. 1 1 . — Comparison of the median mass-size relation found in the PM2GC (continuous black lines), in WINGS (red dotted lines) and those given by Shen 
et al. (2003) for the Sloan r band (green dashed lines). The comparison is done using: all galaxies (top left panel), early-type (T < 0) PM2GC and WINGS 
galaxies with n > 2.5 Shen's galaxies (middle left), late-type (T > 0) PM2GC and WINGS galaxies with n < 2.5 Shen's galaxies (middle right), n > 2.5 PM2GC, 
WINGS and Shen's (bottom left), and n < 2.5 PM2GC, WINGS and Shen's (bottom right). The five lines for each dataset represent the median (middle line), \a 
and 2cr relations. All WINGS radii in this figure are based on B-band imaging. Shen et al. (2003) do not provide a mass-size relation for all galaxies. 



older than 8 or 10 Gyr (red). 

We conclude that both in the field and in clusters, as one 
considers galaxies with older LW ages, the median mass-size 
relation shifts to smaller sizes, but this effect is more pro- 
nounced in clusters than in the field. On average, in the mass 
range 3 x 10 10 -2 x 10 n M©, the sizes of "old" (8+ Gyr in 
LW age) galaxies are a factor ~ 1 .2 in the field and ~ 1 .5 in 
clusters smaller than those of "young" (4Gyr) galaxies. 

This effect needs to be taken into account when compar- 
ing the sizes of high-z and low-z galaxies to infer an evolu- 
tion in the sizes of individual galaxies. Our results show that 
selecting already passive galaxies at high-z means selecting 
the most compact ones. Assuming they do not resume the 
star formation activity at later times, the high-z sizes should 
be compared with the sizes of the oldest (in terms of their 
luminosity-weighted age) low-z galaxies, to avoid a signifi- 
cant progenitor bias. 

Moreover, we have shown that compact massive galaxies 
are much more common in clusters than in the field today, 
and that cluster SDGs are older than the field ones. It is there- 
fore logical to expect that a large number of the high-z mas- 
sive galaxies selected to be passive are found in clusters today. 
This is supported also by the expectations obtained from hier- 
archical simulations, which we now discuss. 

5.1. Theoretical expectations 



The median galaxy mass of the ensemble of high-z datasets 
used in this paper is M t = 10 10 % M Q . In simulations, the most 
massive galaxies at high redshifts are located in the extreme 
peaks of the density distributions, and it is reasonable to ex- 
pect that many of them will end up in the most massive struc- 
tures today. 

We use the Millennium Simulation (Springel et al. 2005), 
populated using the semi-analytic model presented in De Lu- 
cia & Blaizot (2007). These simulations show that 40% 
(28%) of all galaxies more massive than M* = 10 n M Q (M* = 
5 x 10 10 M Q ) at z=2-2.5 end up in haloes more massive than 
1O 14 M (WINGS-like clusters) by redshift 0. This fraction 
decreases to 30% (23%) and 20% (16%) for galaxies with 
> 1O U M (M* = 5 x 1O 1O M ) at z = 1.6 and z = 0.6, re- 
spectively. 

Interestingly, if one considers galaxies with M* > 10 n M Q 
(M t = 5x 10 10 M Q ) and already passive (SFR < lM Q /yr, 
all with a specific star formation rate SSFR = SFR/M+ < 
lO -11 }^" 1 ) at z=2.2, 1.6 and 0.6, the fraction that will be in 
> 10 14 M Q clusters at z=0 rises to 59% (45%), 44% (36%) 
and 25% (22%), respectively. 

Hence, a large fraction of high-z massive galaxies, and the 
majority of high-z passive and massive galaxies have evolved 
into today's cluster galaxies. Clusters at low redshift are the 
most natural location to look for the descendants of the major- 
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FIG. 12. — Median mass-size relation for PM2GC (left) and WINGS (right) galaxies with LW ages > 8 Gyr (red continuous lines) and LW ages < 4 Gyr (blue 
continuous lines). Upper and lower quartiles are shown as errorbars. For WINGS, dotted lines refer to LW ages < 6 Gys (blue) and > 10 Gyr (red). 



ity of high-redshift massive and passive (and compact) galax- 
ies. 

From simulations it is also possible to obtain a prediction 
for the fraction of high-z massive and passive galaxies that 
will end up being the central galaxy of today's clusters, there- 
fore the BCG of clusters like WINGS's. 

Cluster BCGs are the only type of galaxies for which a large 
evolution in size was required by our WINGS study (V10, see 
also Bernardi 2009). The comparison of WINGS clusters and 
EDisCS clusters showed that the mean size and mass of BCGs 
have respectively increased by factors of ~ 4 and ~ 2 between 
Z ~ 0.7 and z ~ 0.04 (Valentinuzzi et al. 2010b). 

Simulations predict that 40% (23%), 26% (14%) and 7% 
(3%) of galaxies more massive than logM+ = lO n M (M t = 
5 x 1O 1O M ) and passive (SFR < lM Q /yr) at z=2.2, 1.6 and 
0.6 become BCGs of clusters at z=0. Hence, a significant frac- 
tion of massive and passive (and compact) galaxies at z > 2 



should be the progenitors of cluster BCGs. As they are sitting 
in the best location of the Universe for accreting other galax- 
ies, their expected minor merging rate is extremely high, and 
the time between z = 2 and z = could be sufficient for a colos- 
sal size growth to reach the > 20 kpc sizes they have today. 

5.2. Progenitors and descendants: the amount of size 
evolution in massive galaxies 

A summary of our results regarding the relations between 
mass, size, age and environment is presented in Fig. [13] The 
most important trends for the LW age are exemplified in the 
left panel: 

a) at a given size, more massive galaxies have older LW 
ages, both in clusters and in the field. The mass of a galaxy 
has a strong influence on the epoch when it stops forming 
stars. 

b) at a given mass, galaxies with smaller radii have older 
LW ages, both in clusters and in the field. The size (compact- 
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ness) of a galaxy influences the epoch when it stops forming 
stars. 

c) at given mass and size, galaxies in clusters have older 
LW ages than galaxies in the field. On top of the dependence 
on mass and size, environment has a significant effect on the 
epoch when galaxies stop forming stars. 

When trying to quantify the impact of stellar-age selection 
effects on the apparent evolution of the mass-size relation of 
passive galaxies, the "age" to be considered is the luminosity- 
weighted age, because the spectrophotometric properties that 
are used by the high-z studies to select "passive" galaxies 
(color, absence of emission lines, SED age etc) are related to 
the LW age, that reflects the epoch of the last star formation 
episode. It is however interesting to analyse also the relations 
between mass, size, age and environment considering the MW 
age, related to the epoch when the bulk of stars were formed. 
As shown in Fig. [13] 

d) at a given size, more massive galaxies have older MW 
ages, both in clusters and in the field, though the differences 
of age with mass are smaller than for the LW age, at the level 
of 2 Gyr at most in our mass range. The mass of a galaxy has 
a significant influence on the epoch when galaxies form the 
bulk of their stars. 

e) the trends of MW age with radius at fixed mass are gener- 
ally weaker than for the LW age, especially for the two lowest 
mass bins. The size of a galaxy seems to have a rather weak 
influence on when the bulk of stars were formed. Our capa- 
bility to detect such influence is of course hampered here by 
the fact that the age resolution of spectrophotometric models 
at such old ages is poor. 

f) there is a noticeable dependence of MW age on environ- 
ment, with cluster galaxies of a given mass and size being 
older than their field counterparts. On top of the dependence 
on galaxy mass and size, environment plays a significant role 
in determining when galaxies form the bulk of their stars. 

The results we have presented highlight the necessity to 
take into account the dependence of the mass-size relation 
on galaxy luminosity-weighted ages and environment when 
studying the evolution of galaxy sizes. Failing to do so has 
the effect of overestimating the evolution of the sizes of indi- 
vidual galaxies. 

To quantify by how much the size evolution could be over- 
estimated, we compare the high-z median mass-size relation 
(solid black line in Fig. [14]! with that of old (LW age> 10 
Gyr) WINGS and PM2GC galaxies (red lines in figure). We 
find an average evolution of 0.2-0.25dex, i.e. a factor 1.6-1.8, 
for WINGS and PM2GC respectively!!] This factor represent 
by how much the median mass-size relation shifts between 
high-z and low-z when considering the most likely descen- 
dants of the high-z progenitors. We do not attempt to assess 
the amount of evolution as a function of galaxy mass (see e.g. 
Ryan et al. 2012 and Huertas-Company et al. 2012), because 
this would require a mass-representative high-z sample. 

Comparing instead the high-z relation with the PM2GC re- 
lation for all galaxies with n > 2.5 or all of today's passive 
galaxies at face value, one would find an evolution of about 
0.5dex, a factor 3-3.2. Therefore, environmental and selection 
effects can account for half of the apparent effect. 

We conclude that a significant fraction of the massive com- 
pact galaxies found at high-z has evolved into today's cluster 

9 We note that this amount of evolution is in agreement with the factor 
1.5 at most estimated by V10 comparing the sizes of galaxies of different 
luminosity-weighted ages with galaxies at different redshifts. 



compact galaxies. When accounting for the age-selection ef- 
fects and the hierarchical assembly bias, the evolution of the 
high-z population in size is smaller than would be inferred at 
face value considering as descendants the whole passive low-z 
population regardless of age and environment. The size evo- 
lution of individual galaxies is modest, of the order of a factor 
1.6- 1.8 on average. 

6. COMPARISON WITH PREVIOUS WORKS 

Most of the high-z works use as local comparison the me- 
dian mass-size relation derived for n > 2.5 galaxies by Shen 
et al. (2003) from Sloan r-band imaging, finding several high- 
z points to lie at more than 2cr from such relation. In §5 we 
have seen that our field relation for n > 2.5 galaxies is shifted 
by about lcr to larger radii compared to the Shen et al. rela- 
tion, and that such offset might result from a combination of 
effects (bluer band used in this work, systematic error in size 
estimates used by Shen, etc). 

A few works have tried to identify massive compact galax- 
ies, similar to the high-z galaxies, at low redshifts. Two 
works in particular have searched for superdense galaxies in 
the nearby Universe, both of them using Sloan-based size and 
mass measurements. 

Trujillo et al. (2009) searched for galaxies with M* > 
8 x 10 10 and R e < 1.5 kpc in the Sloan at < z < 0.2 us- 
ing NYUVAC radii and mass measurements. After exclud- 
ing stars, QSOs and cleaning the sample from those close to 
bright stars, in pairs and edge-on disks, they were left with 29 
SDGs. Applying Trujillo et al. (2009) criteria to define SDGs, 
we find 4 galaxies in our sample, out of 295 ( 1 .4%) more mas- 
sive than Trujillo's limit. As we find in this work and as V10 
found in nearby clusters, Trujillo's SDGs are mostly fast ro- 
tators SOs (Ferre-Mateu et al. 2012, Trujillo et al. 2012). 
The most intriguing result of their work is the fact that their 
SDG stellar population ages are young, suggesting that these 
galaxies have experienced recently enormous bursts of star 
formation (Trujillo et al. 2009, Ferre-Mateu et al. 2012). In 
contrast, the stellar populations of our SDGs selected accord- 
ing to Trujillo's criteria are old: the MW age is always greater 
than 8 Gyr, the mean LW age is 6.3 Gyr and is > 5 Gyr for all 
except one galaxy with LW age=1.2Gyr. 

Another work searching for SDGs in the local Universe 
selected SDSS red sequence galaxies with > 10 ia7 M Q 
with sizes less than half the size predicted by the Shen et al. 
(2003) n > 2.5 relation, excluding galaxies that were disk- 
like by visual inspection, showed marked asymmetries, were 
merging or affected by bright stars (Taylor et al. 2010). We 
note that excluding visually disk-like galaxies they probably 
would have excluded most of our sample. Using de Vau- 
coulers DR7 radii and MPIA SDSS stellar masses, they found 
63 galaxies at 0.066 < z < 0.12. Using V10 criteria to de- 
fine SDGs, Taylor et al. get an SDG fraction of 1.3% for 
red sequence galaxies. If we use Taylor et al. (2010) crite- 
ria to define SDGs (red sequence galaxies, M t > 1O 1O 7 M , 
log(R e )(kpc) < 0.56 *{logM i< - 9.84) -0.3), but without ex- 
cluding disks, we identify 13 galaxies out of 589 above this 
mass limit, therefore a fraction equal to 2.2%. Excluding SOs 
and spirals, only 3 galaxies remain, a fraction 0.5%. In con- 
trast with Trujillo's findings and in agreement with ours, Tay- 
lor's SDGs have old stellar population, with LW ages typi- 
cally in the range 6-10 Gyr. 

Overall, our results agree with previous works in finding 
that superdense galaxies are a small fraction of the global 
general field population at low redshift, the exact value of 
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FIG. 13. — Median luminosity-weighted (left) and mass-weighted (right) age of galaxies as a function of effective radius for four bins of galaxy masses. Solid 
lines: field (PM2GC). Dashed lines: clusters (WINGS). The rms of the median is shown as errorbars only for the lowest mass bin, for clarity. Only points with 
more than 3 galaxies are plotted. 



such fraction depending of course on the adopted definition 
of what is superdense. It remains however to be understood 
why the properties of galaxies selected applying Trujillo's cri- 
teria on the SDSS are so much different from those found with 
Taylor's and our selection (see also discussion in Taylor et al. 
2010 about the comparison with Trujillo et al. 2009). 

In a separate work, Trujillo et al. (2011) investigated the 
issue of the age dependence on galaxy size using a sample 
of visually classified ellipticals from the SDSS, and an HST 
GOODS sample at higher-z. They found no trend of size with 
MW age, at fixed mass, both in the SDSS and in the z ~ 1 sam- 
ple. This agrees with the consistency between the MW age of 
SDGs and all galaxies that we also find in Fig. [7] However, 
these authors also mention that their result does not change if 
they used LW ages. Interestingly, they do find a variation of 
size with MW age at fixed dynamical mass, in the direction of 
older galaxies being more compact, as we find. 

It should be noted that they selected visually morphologi- 
cally classified ellipticals (no SOs included), while the bulk of 
the age selection effects should show up when including pas- 
sive galaxies with a disk component, especially SOs that are 
a very important component of the passive/early-type popu- 
lation, representing ~ 50% and 60% of the total early-type 
massive galaxy population in the general field and in clus- 
ters at low-z, respectively (Calvi et al. 2012). Disky galax- 
ies are common among high-z superdense massive galaxies 
(van Dokkum et al. 2008, Cimatti et al. 2008, van der Wei 
2011, and other high-z references - see also Chevance et al. 
2012 for the possibility that structurally high-z massive com- 



pact galaxies are unlike those of any local galaxy sample), and 
they dominate the SDG populations at low-z both in the field 
and in clusters, as shown in Table 2. It is also important to 
stress that the size-age trend at fixed mass is much more pro- 
nounced in clusters than in the general field, as we discussed 
in §5, therefore it is harder to detect in the field than in clus- 
ters. 

Finally there are a few works in the literature investigat- 
ing the environmental dependence of the relation between 
galaxy size and mass, none of them in the local Universe. The 
lowest redshift study is the one of Maltby et al. (2010) at 
Z ~ 0.17, who investigated the stellar-mass-size relation for 
ellipticals, lenticulars and spirals, separately, in the A901/2 
STAGES field. Using major axis effective radii and COMBO- 
17 masses, they analyzed a cluster, a cluster core and a field 
sample, the latter selected in the cluster area on the sky in red- 
shift intervals that avoid the cluster. For ellipticals, lenticulars 
and high-mass spirals they found no evidence for an environ- 
mental dependence, while for intermediate/low mass spirals, 
field galaxies are larger at 2a than in the cluster. An analysis 
of the mass-size relation for each morphological type is be- 
yond the scope of this paper, but in future works it would be 
interesting to assess how the environmental dependence of the 
size-mass relation is related to the variation of morphological 
mix with environment. 

Going to higher redshifts, apparently contrasting results 
have been found. Raichoor et al. (2012) found that visually 
classified early-type galaxies, mostly atM* < 2 x 1O U M0, are 
more compact in the Lynx supercluster at z = 1.3 than in the 
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FIG. 14. — Median mass-size relation at high-z and low-z. Red lines: WINGS 
points. Green lines: PM2GC all galaxies, with 1 and 2cr. Large symbols indicate 

field at similar redshifts, with a 35-45% of the cluster pop- 
ulation being SDGs according to our definition. The evolu- 
tion observed for the cluster mass-size relation is consistent 
with the hypothesis that larger size, bulge-dominated spirals 
transform into earlier types and provoke the shift of the early- 
type mass-size relation with redshift (Mei et al. 2012). In 
contrast, Rettura et al. (2010) found comparable mass ver- 
sus size relations in red sequence galaxies in a cluster and 
in the field at z ~ 1.2. Moreover, quiescent (by a combina- 
tion of colors) galaxies in a proto-cluster at z ~ 1 .6 have been 
found to have larger average effective sizes compared to qui- 
escent field galaxies at fixed mass, resulting in a deficit of 
compact galaxies and suggesting an accelerated evolution in 
clusters compared to the field (Papovich et al. 2012). Even 
more recently, Huertas-Company et al. (2013) do not detect 
an environmental dependence of the galaxy mass-size rela- 
tion and the size growth in the field-to-group halo mass range 
(< 10 14 M Q ) in COSMOS and SDSS. The environmental de- 
pendence of the mass-size relation at high-z is clearly still an 
open issue. 

7. SUMMARY 

In this work we have analyzed the Padova Millennium 
Galaxy and Group Catalogue to search for SDGs in the lo- 
cal Universe, study their properties and investigate how their 
fraction and characteristics change with environment, using 
as comparison the WINGS sample of cluster galaxies. Our 
main results can be summarized as follows. 

SDG galaxies, defined to have densities > 



and PM2GC galaxies with LW ages > 10 Gyr. Black line: median to all high-z 
high-z datapoints. 

3 x 10 9 M Q kpc" 2 , represent 4.4% of all galaxies at masses 
3 x 10 1() M© < M t < 4 x 10 U M©. This corresponds to a 
number density of 4.3 x 10~ 4 h 3 Mpc" 3 . If we used literature 
radii and mass estimates instead of our own, we would 
get similar numbers. The SDG radii and mass densities 
distributions resemble those of high-z massive and compact 
galaxies. 

SDGs tend to be flattened, with an average intermediate 
Sersic index (< n >= 2.8), consistent with the fact that 70% 
are SOs. They are generally red and passive, with a median 
LW and MW ages of 5.4 and 9.2 Gyr, respectively. Their 
velocity dispersions and dynamical mass estimates are con- 
sistent with them being exceptionally dense compared to the 
general galaxy population at similar masses. 

The incidence of SDGs and their luminosity-weighted ages 
depend strongly on environment. The SDG fraction is 3 times 
higher and the median age 4 Gyr older in clusters above 
500km s" 1 than in the general field. Other SDG properties, 
such as morphological distribution and median Sersic index, 
are very similar for cluster and field SDGs. We find no evi- 
dence for a trend with local galaxy density, neither a change 
with global environment going from isolated galaxies to bi- 
nary systems to groups. 

We study how the mass-size relation changes with the stel- 
lar age of the galaxy and with environment. We find that the 
mass-size relation shifts to smaller radii for older luminosity- 
weighted galaxies. This effect is mild in the field, and much 
stronger in clusters. 
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ID 


R e (pixel) 


n 


b/a 


61514 


9.553 


7.797 


0.488 


16076 


13.133 


0.854 


0.678 


55832 


9.457 


3.295 


0.532 


22663 


9.844 


4.566 


0.947 


43503 


8.667 


3.111 


0.779 


18454 


5.676 


3.751 


0.332 


37797 


10.338 


1.242 


0.576 


18855 


5.477 


3.544 


0.426 


12239 


7.964 


3.887 


0.592 


19977 


6.088 


3.251 


0.618 


24151 


5.372 


2.096 


0.413 



TABLE 5 

Surface photometry parameters computed with GASPHOT for the PM2GC mass-limited sample. Columns : ( 1 ) PM2GC ID number 
as in Calvi et al. (201 1); (2) Galaxy maior axis radius in pixel. The pixel size is 0.333 arcsec; (3) Sersic index; (4) Axial ratio b/a. 
The full catalog is available as online material at the iournal and at the CDS. 



Moreover, cluster SDGs with old luminosity-weighted ages 
represent a significant fraction of all old SDGs at low-z. Un- 
fortunately we can only place a hard lower limit to the fraction 
of SDGs that are in clusters, which is 17%. These old SDGs 
have sufficiently old stellar populations to be the descendants 
of the high-z galaxy population assuming they have not re- 
sumed a star formation activity at later times. 

Based on these results and supported by the expectations 
of hierarchical simulations that predict that 60% of galax- 
ies that are already massive (> 10 n M Q ) and passive (SFR 
< IMoyr" 1 ) at z ~ 2 end up in clusters at z = 0, we argue that, 
in order to assess the amount of size evolution of individual 
galaxies, it is more appropriate to compare the high-z data 
with old, cluster SDGs at low-z than with the general field 
population of all ages. 

An important fraction of the high-z galaxies is expected to 
be the progenitors of today's cluster BCGs, which in fact have 
undergone an evolution in size and mass that is much stronger 
than that of all other galaxies. For these galaxies, the evolu- 
tion has been a factor 4 in size and 2 in mass between z ~ 0.7 
and z ~ 0.05 (Valentinuzzi et al. 2010b). These are of course 
the most massive galaxies locally, usually above the mass 
limit we adopt here (M* > 4 x 10 n M Q , see Fig. 5 in Valentin- 
uzzi et al. 2010). This is consistent with a mass-dependent 
size evolution found by other studies (Huertas-Company et 
al. 2012, Ryan et al. 2012), which we speculate might be due 
to BCGs. 

When taking age and environmental effects into account, 
we find that the average evolution in size of individual, non- 
BCG galaxies of masses 3 x 1O 1O M < M, < 4 x 1O H M is 
very mild, a factor ~ 1 .6. 



One of our most intriguing results is the dependence of the 
LW and MW ages on galaxy mass, size and environment. The 
LW age, that can be considered a sort of "quenching age" ap- 
proximately equal to the time elapsed since galaxies stopped 
forming stars, is significantly influenced by mass, size and 
environment. Both in clusters and in the field, at a given size 
more massive galaxies have older LW ages and, at a given 
mass, smaller galaxies have older LW ages. At fixed mass and 
size, cluster galaxies have older LW ages than field galaxies. 
The existence of a link between mass, structure and stellar 
age was already discussed in D'Onofrio et al. (2011) based 
on different arguments. 

The MW age, that represents the "true" average age of stars 
in galaxies and is more related to the epoch when galaxies 
formed the bulk of their stars, depends strongly on galaxy 
mass and environment, and less strongly on galaxy size. 
These findings may help guiding models of galaxy forma- 
tion and evolution to identify the mechanisms determining the 
birth epoch of most stars and the quenching epoch of star for- 
mation in galaxies. 
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APPENDIX 

GALAXY SIZE, AXIAL RATIO AND SERSIC INDEX CATALOG 

We release the catalog of sizes, axial ratios and Sersic indices measured with GASPHOT for the PM2GC mass-limited galaxy 
sample used in this paper. Table Al shows a subset of the catalog, available in full as online material at the journal and at the 
CDS. The catalog of positions, galaxy stellar masses, redshifts, environment and other galaxy properties can be found in Calvi 
et al. (2011). The morphological catalog is released with Calvi et al. (2012, submitted), and upon request to the author before 
publication. 

GIM2D radii, axial ratios and Sersic indices used in this paper are given in Allen et al. (2006) and can be found at the 
Millennium Galaxy Catalogue web page http://www.eso.org/^jliske/mgc/, file gim2d_sersic0.92. Images are also released by 
the MGC, following the instructions at the same web page. 
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